The Solar pp and hep Processes in Effective Field Theory 
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Abstract. The strategy of modern effective field theory is exploited to pin down accurately the flux S factors 
for the pp and hep processes in the Sun. The technique used is to combine the high accuracy established in 
few-nucleon systems of the "standard nuclear physics approach" (SNPA) and the systematic power counting of 
chiral perturbation theory (ChPT) into a consistent effective field theory framework. Using highly accurate wave 
functions obtained in the SNPA and working to N 3 LO in the chiral counting for the current, we make totally 
parameter- free and error-controlled predictions for the pp and hep processes in the Sun. 



In this talk, we report on the result of the program sketched in [2] which was made possible by the collaboration with 
L.E. Marcucci, R. Schiavilla, M. Viviani, A. Kievsky and S. Rosati, recently summarized in [3]. 



INTRODUCTION 

One of the ultimate goals of quantum chromodynamics (QCD) in nuclear physics is to make precise model- 
independent predictions for certain processes that figure importantly in astrophysics and cosmology. The only tech- 
nique presently available to achieve such goals in the low-energy domain is effective field theory (EFT) [1]. In this talk, 
we apply a variant of nuclear EFT developed by the authors to the solar pp and hep processes within the framework 
sketched in [2]. In doing so, we rely on the accurate wave functions constructed by Marcucci et al [4] as the essential 
ingredient of EFT coming from the standard nuclear physics approach (SNPA). The power of the proposed scheme is 
the ability to correlate the beta decay processes of A = 2,3,4 nuclei that allows us to fix one unknown constant in the 
theory, rendering possible a totally parameter-free prediction of the two-nucleon and four-nucleon processes. 
The most abundant source of solar neutrinos (carrying 91 % of the total flux) is the pp process 

p + p^d + e + +v. (1) 

This process has been carefully studied [5, 6, 7, 8], and the calculated transition strength, governed by the leading- 
order Gamow-Teller (GT) operator, is believed to be reasonably reliable. However, given its extremely important role 
in the solar burning process, the pp process invites further elaborate studies. Meanwhile, the hep process 

3 He + p^ 4 He + e++v. (2) 

produces the highest energy solar neutrinos, E™ x (hep) — 20 MeV. While the hep neutrino flux is estimated to be 
much smaller than the 8 B neutrino flux, there can be significant distortion of the 8 B neutrino spectrum at its higher 
end if the hep S-factor is much larger than the existing estimates. This change can influence the interpretation of 
the results of a recent Super-Kamiokande experiment that have raised many important issues concerning the solar 
neutrino problem and neutrino oscillations [9, 10]. To address these issues, a reliable estimate of the hep cross 
section is indispensable. Its accurate evaluation, however, has been a long-standing challenge for nuclear and hadron 
physics [11]. The difficulty involved is reflected in the pronounced variance in the documented estimates of the hep S- 
factor. For example, the first estimate given by Salpeter [12] was S(hep) = 630 x 10~ 20 keV-b, which was eventually 



replaced by the much smaller (so-called "standard") value, ~ 2 x 10 keV-b [13, 14]; the latest, most elaborate 
estimation gives 9.64 x 10~ 20 keV-b [4]. The reason for the difficulty in making a precise estimation of the hep S- 
factor is multifold. First, the one-body (IB) GT matrix element for the hep process is strongly suppressed due to the 
symmetry properties of the orbital wave functions of the initial and final states. The main orbital wave function for 4 He 
has [4] symmetry (totally symmetric) under the particle exchanges, whereas the dominant 3 He +p orbital wavefunction 
has [3 1 ] symmetry; the [4] component is forbidden by the Pauli principle when there are three protons. Then, the main 
components of the initial and final states - with different symmetry properties in orbital space - cannot be connected 
by the leading-order (LO) GT operator, which does not contain orbital variables. This means that the non-vanishing 
IB GT matrix element (for the hep process) is due to either minor components of the wavefunctions or higher order 
corrections of the GT operators. Since these are all quite small, the IB matrix element becomes comparable to multi- 
body corrections, e.g., meson-exchange-current (MEC) contributions. A further complicating feature is that there is 
a substantial cancellation between the IB and two-body (MEC) contributions, which can amplify the errors. Finally, 
with the "chiral-filter mechanism" rendered ineffective, at non-vanishing leading order of effective field theory (see 
below), the many-body corrections contain short-ranged operators the strengths of which are not known a priori and 
hence difficult to control. 

The objective of our present work is to make accurate effective field theory (EFT) predictions on the pp and hep 
processes within a single framework. For this purpose, we adopt the strategy that exploits the power of both SNPA and 
heavy-baryon chiral perturbation theory (HBChPT), which is a well-studied low-energy EFT. In this approach, EFT 
enters in the calculation of relevant transition operators. We will calculate them up to next-to-next-to-next-to-leading 
order (N 3 LO); all the operators that appear up to N 3 LO will be considered. To obtain the corresponding nuclear matrix 
elements, we need highly accurate nuclear wave functions. Although it is- at least in principle - possible to also derive 
nuclear wave functions to the appropriate order from HBChPT, we choose not to do so. Instead, we use realistic wave 
functions obtained in the standard nuclear physics approach (SNPA). The potentials that generate such wave functions 
are supposed to contain high orders in the chiral counting, presumably much higher than what can be accounted for 
in the irreducible vertex for the current operators. For a review of SNPA, see Ref. [15]. Such an approach- which is 
close in spirit to Weinberg's original scheme [16] based on the chiral expansion of "irreducible terms"-has been found 
to have an amazing predictive power for the n + p — ► d + y process [17, 18, 19]. 

The basic advantage of EFT is that the SNPA and HBChPT can be combined into a model-independent framework 
based on the first principle. A systematic expansion scheme of EFT reveals that, for the GT transition for which the 
"chiral-filter mechanism" [2] is rendered inoperative, the one-pion-exchange (OPE) operators and the leading short- 
ranged operators have the same chiral order, and hence their contributions should be comparable. This nullifies the 
intuitive argument that the long-ranged OPE contribution should dominate the MEC corrections. The existing SNPA 
calculation, however, lacks this leading short-ranged operator while containing instead some higher order short-ranged 
contributions. What our EFT manages to do is to account for this short-ranged contribution in a way consistent with 
renormalization group invariance. This is a novel way of understanding the so-called "short-range correlation" in 
SNPA. We will see that this aspect indeed plays an essential role both for pp and hep processes. 

Briefly, our approach to HBChPT is as follows. We take only pions and nucleons as pertinent degrees of freedom. 
All others have been integrated out, and their dynamical roles are embedded in the higher-order operators. In the 
scheme relevant to us, it suffices to focus on "irreducible graphs" according to Weinberg's classification [16]. Graphs 
are classified by the chiral power index v given by v = 2 (A — C) + 2L + V„ where A is the number of nucleons 
involved in the process, C the number of disconnected parts, and L the number of loops. The chiral index, v,, of the z'-th 
vertex is given by V,- = di + e, + n,/2 — 2, where <f e, and n, are respectively the numbers of derivatives, external fields 
and nucleon lines belonging to the vertex. The Feynman diagrams with a chiral index v are suppressed by (g/A z ) v 
compared with the leading-order one-body GT operator, with Q standing for the typical three-momentum scale and/or 
the pion mass, and A x <~ ~ 4nf n is the chiral scale. The physical amplitude is then expanded with respect to v. 



THEORY: GT OPERATORS UP TO WLO 

The LO and N 2 LO contributions come from the well-known one-body currents; the GT operator for the a-th isospin 
component reads 



MB 



8 a L T 



1=1,2 



PiOrPi-QiPi 

2m 2 N 



(3) 



with p, = (p/ +p/)/2 and gA ~ 1.267. Corrections to the above IB operators are due to MECs, which start at N 3 LO. 
In our work we include all the contributions up to N 3 LO. We emphasize in particular that, up to N 3 LO, only two-body 
(2B) currents enter, three-body currents appearing only from N 4 LO; thus there is no arbitrary truncation involved here. 
The N 3 LO 2B currents consist of the OPE and contact-term (CT) parts, 



A^b=A^b(OPE)+A^ b (CT). (4) 



The OPE part is given as 
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where % a Q = (x\ xif, with 



x,—, and similarly for o . Since the couplings C34 are determined from the %N 
data [20], £3 = —3.66 ±0.08 and £4 = 2.11 ±0.08, there is no unknown parameter here. Now the CT part is given as 



A!b(CT) = --^2 [Ji(Tfai +^o 2 ) + i2T fl x Ox] , 
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which contains two low-energy constants. However, the Pauli principle (or the "L + S + T = odd" rule) combined 
with the fact that the CT term is effective only for s-wave (L = 0) implies that we need only work with one unknown 
constant, d R , defined by 



d R = di+2d2 + -c 3 + -c 4 + 2 ■ 
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1 



1 



(7) 



Furthermore, the same combination also enters into tritium P-decay, /j-capture on deuteron, and v-d scattering. 
Although d R is in principle calculable from QCD for a given scale A , this calculation is not available at present; 
we therefore need to fix d R empirically. Here we choose to determine d R by fitting the tritium p-decay rate, Tp, which 
is accurately known experimentally [7]. Once d R is fixed, our calculation involves no unknown parameters. 

We calculate the matrix elements of the transition operators with state-of-the-art realistic nuclear wave functions for 
A = 2, 3, 4. We employ the correlated-hyperspherical-harmonics (CHH) wave functions, obtained with the Argonne 
vis (Avl8) potential (supplemented with the Urbana-IX three-nucleon potential for the A > 3 nuclei). To control short- 
range physics in a consistent manner, we apply the same regularization method to all the nuclear systems in question. 
Specifically, in performing Fourier transformation to derive the r-space representation of transition operators, we use 
the Gaussian regularization. This is equivalent to replacing the delta and Yukawa functions with the regularized ones, 



where the cut-off function S\(q 2 ) is defined as 

^ 2 ) = -p(-^). 
The cutoff parameter A characterizes the energy-momentum scale of our EFT. 



(8) 



RESULTS 

The value of d R determined from the experimental value of Tp is 

d R = (1.00±0.07, 1.78 ±0.08, 3.90±0.10) 



(9) 



for the choice of A = (500, 600, 800) MeV, respectively. We list in Table 1 the GT matrix elements for the pp and 
hep processes (in arbitrary units) as a function of A. 



TABLE 1. GT matrix element for the pp and hep processes, calculated for representa- 
tive values of A. The 2B contribution is the sum of the OPE part (^-independent) and the 
CT part (linear in d R ), for each case. 



A (MeV) 


(1B) W 




(2B}« 












500 


4.82 


0.076 


-0.035 d R 


~ 0.041 


-0.81 


0.93 


- 0.44 d R ~ 


0.49 


600 


4.82 


0.097 


-0.031 d R 


~ 0.042 


-0.81 


1.22 


-0.39 d R ~ 


0.52 


800 


4.82 


0.129 


- 0.022 d R 


~ 0.042 


-0.81 


1.66 


-0.27 


0.59 



The results 

We observe that, while the OPE part by itself has a sizable A-dependence, the net amplitude is completely A- 
independent. In other words, the A-dependence of the OPE part has been perfectly removed by that of CT part. The 
A-independence of the physical quantity, (IB) + (2B), which is in conformity with the general tenet of EFT, is a crucial 
feature of the result in our present study. 

The relative strength of the two-body contribution as compared with the one-body contribution is 

5fg ee (2B) PP /<1B) W = (0.86 ±0.05) %. (10) 
This ratio is consistent with the latest SNPA calculation [7], 8^ = (0.5 ~ 0.8) %. The resulting pp 5-factor is 

S(pp) = 3.94 (1 ±0.15 % ±0.1 %) x 10~ 25 MeV-barn, (11) 
where the first and the second uncertainties come from one- and two-body contributions, respectively. 



The hep results 

The general tendency here is quite similar to the pp case; the variation of the two-body GT amplitude is only ~10 
% for the entire range of A under study. The A-dependence in the total GT amplitude becomes more pronounced by a 
strong cancellation between the IB and 2B terms, but this amplified A-dependence still lies within acceptable levels. 

TABLE 2. Contributions to the 5-factor (in 
10~ 20 keV-b) from individual initial channels cal- 
culated as functions of A. The last column gives 
the results obtained in MSVKRB. 



A (MeV) 


500 


600 


800 


MSVKRB 


% 


0.02 


0.02 


0.02 


0.02 




7.00 


6.37 


4.30 


6.38 


3 Po 


0.67 


0.66 


0.66 


0.82 


l Pi 


0.85 


0.88 


0.91 


1.00 


'Pi 


0.34 


0.34 


0.34 


0.30 


'Pi 


1.06 


1.06 


1.06 


0.97 


Total 


9.95 


9.37 


7.32 


9.64 



Table 2 shows the contribution to the 5-factor from each initial channel, at zero cm. energy. For comparison we 
have also listed the latest results based on SNPA [4] (which we refer to as MSVKRB). It is noteworthy that for all the 
channels other than 3 Su the A-dependence is very small (< 2 %). While the GT terms are dominant, the contribution 
of the axial-charge term in the 3 5i channel is sizable even though it is kinematically suppressed by the factor q. 

The results given in Table 2 lead to a much improved estimate of the hep 5-factor: 

S(hep) = (8.6 ±1.3) x 10- 20 keV-b, (12) 

where the "error" spans the range of the A-dependence for A = 500-800 MeV. This result is to be compared to that 
obtained by MSVKRB [4], S = 9.64 x 10~ 20 keV-b. To decrease the uncertainty in Eq.(12), we need to reduce the A- 
dependence in the two-body GT term. According to a general tenet of EFT, the A-dependence should diminish when 
higher order terms are included. A preliminary study indicates that it is indeed possible to reduce the A-dependence 
significantly by including N 4 LO corrections. 



DISCUSSION 



By determining the only parameter of the theory d from the experimental data on triton beta decay, we have succeeded 
in making rather accurate EFT predictions (up to N 3 LO) in a parameter-free manner for both pp and hep processes. 
These predictions turn out to give support to the latest SNPA results. 

The prediction for the pp prediction comes out to be independent of the cutoff scale A, which means that it is fully 
consistent with the tenet of EFT. On the other hand, there remains some A-dependence for the hep process, which 
could be due to many-body nature absent in the pp case. Even so, it is remarkable that the theoretical uncertainty can 
be reduced from "orders of magnitude" to ~ 20 %. This uncertainty can be further reduced if N 4 LO terms - which 
involve no additional unknowns - are taken into account. 

We should note that by using the "exact" wave functions, we are sacrificing the strict adherence to chiral order 
counting in favor of predictivity. The counting error committed therein comes at one order higher than that accounted 
for in the irreducible vertex for the current , i.e., at N 4 LO in the present calculation, and should be small for the whole 
scheme to make sense. This can be checked by an N 4 LO calculation for which the counting error would come at N 5 LO 
or higher. Furthermore the notion implicit in our approach that the IB matrix element calculated with the "exact" wave 
functions should be taken as "empirical" could be checked by looking at the hen process 

3 He + «^ 4 He + y 

in which the IB matrix element is expected to suffer the same suppression due to the symmetry properties of the initial 
and final states. 

It seems likely that an EFT calculation of the hep process that adheres to the strict power counting - such as the 
pionless theory [21]- would be obstructed by a plethora of unknown parameters that are difficult to completely control. 
If such a calculation is feasible, however, it would be interesting to see if and how the symmetry suppression and the 
sensitive cancellation encountered in our version of EFT - where an accurate IB matrix element plays a key role - 
could either be circumvented or manifest themselves in the description. 
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